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ABSTRACT 


The purpose of the present note Is to treat some problems of attitude 
stability of spinning satellites in a rigorous manner and to show that, with 
certain restrictions, linearized stability analysis correctly predicts the 
attitude stability of spinning satellites, even In the critical cases of 
Liapunov stability theory^ 
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INTRODUCTION 

Since the attitude instability experience by Explorer 1, many papers 
have been written on the effects of Internal dissipation on the attitude 
stability of spinning satellites. In the engineering literature, stability 
analysis is restricted to the variational or linearized perturbational equations, 
despite the fact that spinning satellites are almost always critical cases 
in the Llapunov-Poincare stability theory. This is certainly true in the case 
of dual spin satellites, which have the further complication that the linearized 
perturbational equations have periodic coefficients. 

The purpose of this note is to treat some problems of attitude stability 
of spinning satellites in a rigorous manner and to show that, with certain 
restrictions, the linearized stability analysis correctly predicts the attitude 
stability of spinning satellites. 


m 
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1 . PetumbllnR of a Spacecraft Using Passive Torsional Dampers 
Formulation of Problem 

Consider a spacecraft which is designed to spin about axis 1, the axis 
of maximum moment of inertia > to provide an artificial gravity field for the 
crew. Attached to the spacecraft on axes 2 and 3 are torsional dampers, conslst~ 


ing of inertia wheels of polar moment of Inertia J^, (i»2,3) with torsional 
springs with restoring torque K^f(B^) and damping torque Let 

l 2 » Ij be the moments of Inertia of the spacecraft about the 1,2 and 3 axes 
respectively, including the moments of Inertia of the dampers 


Let 1«2,3 

Suppose that owing to collision with another spacecraft, which Is at- 
tempting to dock with the first spacecraft, a tumbling motion results. Let 
o)i, W 2 » Wg be the angular velocities of the tumbling motion with respect to 
the body-fixed axes 1,2 and 3 respectively. For the safety and comfort of 
the crew, and to make docking possible, the spacecraft must be detumbled and 
returned to a state of simple spin about the 1 axis. 


Equations of Motion 

• % 

12^2'^ “l^3^^l”^3^'^‘^2®2”'^3^1®3 

I3W3 + WiW2(I2“Ijl)+‘J3*^3‘^J2‘^102 “ ^ 

J2<S2+ W2>+°2®2'^^2^^®2^ =0 

J3(e3+ W3)+D303+K2f(e3) =0 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 

(1.4) 


(1.5) 


( 3 ) 






n 


Where 5 


f(e>— f(- 0 ) 

0f(e)>O Oj^O 


Llm -1 , 

e->o ® 


0 


f(y)dy-F<0)>O OjtO 


( 1 . 6 ) 


!)j,K.>0 1-2,3 

11 


Global Stability 
Let 




) (1.7) 


the function V is clearly positive definite 

V*IiWiWi+1 2‘^2*^2’*'^3‘^3***3’*'‘^2 ^^2^2 '’^2® 2^® 2® 2^^*^2^ ^ ® 2^ ® 2 

+J2(w303 ■^‘^'3®3'^®3®3^"^^3^(®3^ ®3 (1.8) 

Using equations (1.1) , (1. 2) , (1, 3) , (1.4) and (1.5) to evaluate V along 
the trajectories of the motion, we have; 

V«=-[D202+D303]^O (1.9) 

The function V is positive definite and its time derivative along the tra- 
jectories of the motion is negative seml-def inlte, therefore V is a Liapunov 

• 

function and the tumbling motion is globally Liapunov stable. We note that V 

• • 

is only semi-definite and vanishes when e 2 ’= 03 ='O. Equations (1.4), (1.5) show 

f* *i 

that 02 and 0^ are not zero unless. 


<A) 


or 


(a) W2"«j"0 and 02"03"0 

(b) K2f (02)“'"'J2‘*^2 


Examination o£ equations (1.1), (1.2) and (1.3) shows that 

t • • •• • 

bo satisfied In general unless Wj^»a) 2 “W 2**0 , » 

following conditions hold 

I ) , « 2 *‘^ 3“0 

II) w^»W2«0 

III) 


condition (b) cannot 
1"2,3 and one of the 


This set of conditions are simply the equilibrium solutions of the set of 
equations (1.1) ,(1.2) , (1,3) , (1.4) and (1,5). With this exception, 
only on a set of measure zero. Thus, using (1.9), 


V(t)-V(0)»- 


[D202+D3Q3ldt<O 


( 1 . 10 ) 


Hence, the function V(t) decreases along the trajectories of the motion. 

V(t) must therefore tend to a limit corresponding to one of the equilibrium 
solutions. The particular limit to which all motions ultimately tend for large 
time Is determined by the stability of the equilibrium solutions. Clearly all 
motions will tend in the limit to the largest invariant set, which corresponds 
to a stable equilibrium solution. 

Stability of the Equilibrium Solutions 

Examination of equations (1.1) through (1.5) shows that there are three 


equilibrium solutions. 


(5) 


uj^»0 , 1"X,2,3 ; fij*^j**^j**® » J“2,3 

I) w^l*0 , u)2“W3*0 

II) U2f‘0 , Wj^*»W2"0 

III) u^l^O , Wj^*W2*0 

Case (1) 

“l-“lO • “2-V“ 

Let 

“l"“lO+^ V" * “3"' 

Oj-Oj j“2,3 

Perturbing about the steady state solution and retaining only the linear terms 
in the equations of motion, one obtains 


I^5»0 


** ♦ • 

J3(a3+?)+D3i3+K3<x3 ' 



(1.11) 

3“l0“2‘® 

(1.12) 

2“l0“2’‘® 

(1.13) 

«=o 

(1.14) 

bzQ 

(1.15) 


D, 


D 


J. » J„ *'2 


2 -J 

=Pr 

2 


°3 h 2 

J3=^3 ' J 3 -P 3 


Let 



( 6 ) 


Define 


fl 


2 “Io«r’:2><ii-^3) 


’■ 2’'3 


The characteristic equation for the system of linear differential 
equations (1.11) through (1.15) is 



Let 



i“2,3 

If J 2 and Jg are selected such that 



i»2,3 y>0 








(1.16) 


(1.17) 


(1.18) 


<7) 


the characteristic equation (1,16) becomes 

1 {fX^+0X^*fp^A)Wf(l+n)X^+0A+p^]^3 •*0 (1,19) 

which may be writ ten in the form 

fa)-Xgj^(X)B2(X)»0 (1.20) 

where 


g ( X ) - 1 1 X a X ^+p ^ X ] +n [ ( 1+M ) X 6 X+ p ^ ] 
g2(x)-"i[x\(3X^+p^x]-fn[(i+u)x2+3X+p^| 


(1.21) 


Usin)f‘ Ci^iuchy's Principle of the argument, or Nyqulst's criterion, It Is easily 
shown that gj_(X), 1*1,2, have zeros only in the left half X plane. Thus. 

X^»0 , ReX^<0 16(2,7) 

This is clearly one of the criwical cases in Liapunov stability theory, 
however, using Theorem AIII of the appendix, we see that the full perturbation 
equations are Liapunov stable. Thus, the equilibrium solution (1) is stable. 

Case (ii) 

W2*W2o 5^0 , a)^*(D2=0 

Let 

®l““i 


Perturbing about the equilibrium solution and retaining only the linear terms 
in the equations of motion, we have; 


(8) 


Ij^?+aj2QC(l3"l2^'^'^3‘*^20®3''^ (1.22) 

I2n+J2«2 

13^+0)20? (I2“li)«3«i3 -0 (1-24) 

J2(a2+n)+D2“2’^V2 

J 3 ( 613 + 0 + 03(5 3 -»-K 3 a 3 -0 (1.26) 


Let 


2 - 1^13 


D. , K. 

®i“ jT ’ ’’i’ JT 
1 1 




^ (1.27) 




The characteristic equation for this set of linear differential 
equations is: 


p(i- y +e2>+P^]p(i- ^)+63''-e3«^-«2'>3 

+(pr«2(i+i^))^'] ■» 


(1.28) 


Application of Cauchy’s principle of the argument on Nyquist’s stability 
criterion Immediately shows that 


Xl=0 X3<0 

ReX^<0 1=4, 5, 6 , 7 


(1.29) 


since X2>0 , application of Theorem II shows that the full perturbatlonal 
equations are unstable In the sense of Liapunov stability theory. 

%;a8e (111) 

oij-Ujoi'O Uj-Uj-O 

Let 

Perturbing about the equilibrium solution and retaining only the linear terms 
In the equations of motion, we have: 



(1.30) 

l2^“305<^r^3>+-’2“2 

(1.31) 

l3«+^3“3 "0 

(1.32) 

^2 ia2+n)-»>2^2*^2°‘2 

(1.33) 


(1.34) 




> (1.35) 


The characteristic equation for this set of linear differential 


Let 




,2 “ 30^^1 ^ 3^ ^^2 ^ 3 ^ 

2 . 

' 12 


h 2 ^1 

3i= ; p^“ , 1=2,3 


equations is: 


( 10 ) 




(1.36) 
stability 

^ (1.37) 

Since there are two eigenvalues whose real parts are positive, application 
of Theorem All shows that the full perturbatlonal equations are unstable In the 
sense of Liapunov stability theory. Thus we s^e that the only stable equilibrium 
solution Is that corresponding to Case (1) , ( 1 ) 2 * 012 = 0 . From the analysis 

of global stability we know that the function V(t) (1.7) decreases along the 
trajectories of the motion and tends to a limit corresponding to a stable equili- 
brium solution, the only stable equilibrium solution is that corresponding to 
spin about the 1 axis, the axis of maximum moment of Inertia. Thus we have 
shown rigorously that It Is possible to detumble a spacecraft using only passive 
torsional dampers. Edwards and Kaplan (1) have treated the problem of automatic 
detumbling of a spacecraft using the motion of a servo-controlled Internal mass.. 
Their treatment is heuristic rather than rigorous. 

2. Stability of a Dual Spin Satellite 

The stability of dual spin satellites has been examined by a number of 
authors; however, in the case where the rotor and the platform both exhibit 
internal dissipation, the analytical solution was first presented by Sarychev 




Application of Cauchy's principle of the argument, or Nyqulst's 
criterion immediately shows that (1.36) has roots; 


Xj^*0 , ReX^>0 1*2,3 


ReX.<0 1*4, 5, 6, 7 
1 


) 






ill) 


and Sazonov (2) who used Floquet Theory. In this note the effects of Internal 
dissipation will be modelled by torsional dampers in both rotor and platform. 

It will be shown that the linearized stability analysis is rigorously justified 
and it will also be shown that the linearized stability analysis can be performed 
quite simply by using Lagranges method of variation of parameter. 

Formulation of the Problem 

The dual spin satellite consists of two rigid bodies with a common axis 
of rotation (axis 3) 

Let the axial moment of inertia of the rotc/r be J, 

Let the total axial moment of inertia of the satellite be C (rotor plus plat- 
form, plus dampers) 

Let the total equatorial moment of inertia of the satellite be A (including the 
dampers) 

Let Ib and Ib^ be the polar moments of inertia of the damper wheels on the 
platform and rotor respectively 

Let and K 2 be the damping and stiffness parameters of the torsional 

damper on the platform. Let and K 2 ^ be the corresponding parameters for 

the rotor damper: 

Let and ^2 be the angular velocities of the platform with respect to the 

1 and 2 axes respectively. Let be the angular velocity of the platform 

about the 3 axes 

Let ^ be the angular velocity of the rotor about the 3^ axis relative to 
the platform, where the angle (measured about the 3^ axis) defines the 
orientation of the body fixed axes of the rotor With respect to body fixed 
axis of the platform. 


( 12 ) 


let Tg be the frictional torque of the rotor bearings 
Let be the torque of the despin motor. 


Equations of Motion 


J(w2+^)-I|j6 [u)2COs \|i -Wj^sln ij< ]+Tg-Tj^"0 


tp -uj^sln ]-Ij^W24>*0 

Aojj^+[(C-A)w2+Jt|j]w2+l|j[0cos i|) “0(o)2-fi|))sln i|/ 


A(i)2”[ (C-A)wj+Jij;]a)^+I^[0sln ^ ■he(Uj+tJ/)co8 <p 
i|j 0+K^0+K20+I^ { (u)j^+a)2t|^) cos ip ^ ] *0 

I^^J^'+Kj^}+K2^+I^aij^«0 




r (2.1) 


Where 0 are the rotation angles of the torsional dampers on the platform 
and rotor respectively. 


Steady State Solutions 

If the torque of the despin motor just balances the bearing friction 
torque when ip’^a , then the steady state solution is; 

, Jp=a , ip^ot 
0=(j)=0=<f)=a)j^=W2=O 



( 2 . 2 ) 


(13) 


Linearized Stability Equations 
Let 


o 

f T*0t ; r» ~ 

J 0 




0) 


^ ■''l i V "'’2 ' 


''i 

^ A * ^ A ’ *^1 1,0 ’ *^1 I' a 

b b 


K„ K„ 

V s — r — , _ 2 

r% f »>n 


2 ^ 2 * *'2 ^,2 




L (2.3) 


J 


In Appendix C, It Is shown that the steady state solution is stable for 
sufficiently small perturbations, provided the following conditions are satisfied 

3>0 (2.4) 


and the systems of equations: 


T 




(U) 


dc 

dv, 

.1 

dt 


2 r 2 n 

•fXv 2 +W ^2 cos T **(l+r)~ sin tJ “0 

-^V 2 +Mr ^ sin T+(l+r)^ cos ij -0 


[(dT ■^'' 2 ) ■' + (dT" ■'’ 1 ) {] 


^+Kifi+K'j0+f(^+V2l COST + 


^ ‘t’ AV 4.K (j)4. atO 

. 2 dT dT ^ 

dT 


V (2.5) 


is Liapunov asymptotically stable. 

Equations (2.5) have periodic coefficients and may be rewritten in 
standard form as 

d£ 

df “A3<t)E 


( 2 . 6 ) 


where 


and 


AgCT-t tt)=A 2 (t) 


(2.7) 


d^ 

dT 

9 

di 

dt 


( 2 . 8 ) 


(X5) 


The stability of equation (2.6) may then be Investigated by using Floquet 
theory, as was done by Sarychev and Sazonov (1). Alternatively the stability 
of equations (2,5) may be Investigated directly in the case where p are 
small, by using Lagrange's method of variation of parameters. 

If e"Max(p,p^) and e«l , then equations (2.5) are of the type treated 
In Appendix B. 

Let 



cos Xt sin Xt 
sin Xt -cos Xt 


a(T) 


(2.9) 


where 


i(t) = 


'^a^(T) 




Substituting Into equations (2.5) we have: 


da*| dSn Hfi 

cos Xt + sin Xt ^ cos t +p' (1+r)^ sin t -]i — k 

dT dT dT 


da 

dT 


d 3 ^ 

— sin Xt - cos Xt Y sin t (1+r)-^® cos t -yr ^ 


dT^ 


dT 








( 2 . 10 ) 


( 16 ) 


^ ® ^ +K2e-(X-l)[aj^co8(X-l)T-a2Bln(X-l)T] 

[da- da„ 1 

cos(X+1)t+ 8 ln(X+l)Tj 


dt' 


+Kj^ ^ ■fK2<}'"X|aj^Bln Xt 


dT 


•C 


da- da- 

. - cos Xt + sin Xt 
dT dT 


.] 


Using equations ( 2 . 10 ) to solve for da^/dx and da2/dT 


«_y' cos(X-l)T-(l+r)y^ sin(X-l)T 
dT 


-IJ i- 


d(|> 


-y — ^ cos Xt -yr -rr sin Xt 
,2 dT 

dT 


I sin(X-l)T+(l+r)y' ^ cos(X-1)t 

^ dT 


— ^ sin Xt +yr cos Xt 

.2 dT 

dT 


2 2 2 

Clearly, if de/dr , d 0/dT , d0/dT , d <|)/dT are bounded, since e-<<l 


e"Max(y ,y* ) , hence 


dr dT 

N l“2| 


( 2 . 11 ) 


( 2 . 12 ) 


O; 0(e) 


(2.13) 


(17) 


Thu., .j(t) and « 2 (t) are .Xowly varying functions of t. honco In equations 
(2,11) we may neglect the terms da^^/dr , comparison to a^^ and 

We may further treat and a 2 <T) as "constant", provided and Kj^ 

are not too small. 

Thus the "Steady State" solutions of equations (2.11) are; 


0(T)"C(T)cOS(X-l)T+D(T)8ln(X-l)TK)(G) 
■' (t)"E(t)cos Xt +F(T)sln(X-l)T+0(e) 

where 


(2.14) 


C(t)« 


[K2“(X-1)^1^+[kJ^(X-1)]^ 


D(t)» 


E(t)» 


-(X-l)f(X-l)K^a2(T)-(K2-(X-X)^)a^(T)] 

[K2-(X-1)^]^+[K^(X-1)]^ 

-X[(K2-X)^a2(T)+XKj^a^(T)] 

[K2~X^]^+[K^X1^ 


> (2.15) 


F(t) 


-XIXK^a2(T)-(K2-X^)aj,(T)] 

[K2-X^]^+[K^Xj2 


Substituting equations (2.14) into equations (2.12), treating a^(r) and 
a^ix) as constants. Consistent with this, we retain the time averaged coef- 
ficients of C(t) , D(t) , E(t) and F(t) in the resulting equations, thus 


we have 


(18) 


da 


- ■^(<X-I)(X+r)lC+ f[X(X+r)lE40Ce^) 


da 


^^2 . |(X(X-fr)jF+0(?:^) 


(2.16) 


Using equation (2.15) to substitute for C, D, E and V in terms of aj^(*f) and 


a2<x) 


2 

ar > 


dSrt 2 

-jp- -aa2+eaj^+0(e ) 


(2.17) 


where 


a— [X+r] 


■^(X-l)Vj^ 


{K'j^-(X-l)^]^+tK'j^(X-l)f [K2-X^]^+tK^X] 




2 , 2 , 


(2.18) 


3“[X+r] 


[: '3 


i^(X-l)2[K'2-(X-l)^] ^ X^[K2«X^] 


[k'2-(X-1)^]^+[k'j^(X-1)]^ [K2-X'"1''+[Kj^X] 


2 , 2 , 


is simply the matrix of (B-31) . 


(2.19) 


The matrix 

The characteristic equation for the system (2.17) correct to 0(e) is 


(A-a)^+$^-0 


(2.20) 


2 2 2 
A -2aA+a +3 *0 


( 2 . 21 ) 


The condition for stability is that 


a<0 


( 2 . 22 ) 








(19) 


The stability condition may be written 


,<k l ^- 




\^L 

2 , 2 . 


>0 


(2.23) 


where 


a' - -y K^a“l)(A+r) 
A» ^ Kj^A(X+r) 


N 

> (2.2A) 


j 


but 


X. (C-A)r-fJ 


X+r* 


Cr+J 


y (2.25) 


and X-l» 


(C-A)r+(J-A) 




Thus we have the following conditions; 

i) System is asymptotically stable if 

A,a'>0 (2.26) 

li) System unstable if 

AjA'<0 (2.27) 

lii) If A a' ■'0 , stability depends on the quantitative relationship 
between A and A^ 

Iv) In particular, if 0. the spin rate of the platform is zero, i.e. 
r=0 


( 20 ) 


then A* ^ 


(2.28) 


a' , JL y' ( j\ 

^ 2 A 


(2.29) 


v) 


Frequisiitly J<A .*. A^<0 , however by making the dissipation 

In the platform sufficiently large, condition (2.22) can always be 
satisfied 

Provided A>0 , the dissipation In the platform may bo maximized 
2 

by setting K 2 »l ; In this case the condition for stability 
becomes: 


r\ •» * 


/ ft 

^ Kj^(X-l) 'A' 


2 K 


1 [K2-(X-l)^]^+rKj^(X-l)] 


>0 


(2.30) 


In Appendix B It Is shown that the stability treatment presented above 
Is rigorously correct for e«Max(p,p^) sufficiently small. 


Other Problems 

1) The technique above has also been used on the problem treated in 
reference (1) and the results agree exactly. 

2) The technique above has also been applied to the case where the despin 
motor is used In conjunction with the products of inertia terms in the Inertia 
tensor of the platform to obtain stability for the dual spin satellite. 


/ 


( 21 ) 


APPENDIX A 


Definitions 


Given the dynamical system: 
dx 

^ “Ax+f,(x,t) 


x(0)<*£ 

. I If (x,t) 1 1 

Llm ' »0 uniformly In t 

I IeI ko 11*11 


If given any 6>0 there exists an e>0 such that ILcil-^e Implies that 
||x(t)|1^6 , V t>0 , then the trivial solution of A.l is said to be Liapunov 
Stable (L.S.). 


If the trivial solution of A.l is Liapunov stable and in addition 
||x(t)|| tends to zero as t tends to infinity, then the trivial solution of 


A.l is said to be 


:able (L.A.S.). 


)unov ins tat 


If given a <S>0 there exists no G>0 such that ll£il^^ implies that 
|lx(t)||^5 , V t>0 then the trivial solution of A.l is said to be Unstable 
in the sense of Liapunov, 


( 22 ) 


Theorem A1 

If A Is a stability matrix, l.e. if ReA^(A)<0 then A,1 is Liapunov 
asymptotically stable provided that | j^l | is sufficiently small# 

Proof 

Case (i) A noridefective 

There exists a nonsingular matrix T such that T~^AT«A , ReX <0 

1 1 


dz 

^ -Az+£(^,t) 


z(0)« ^T~^c 


£(z,t)-T"^f(Tz,t) 


Urn I I .0 uniformly In t 

I lil ho I Ul I 


V(z)*z*z“ z 


V(z)sz!*'z+z*z 


(2ReA)^+2Re^*_g^(_z , t) 

=-aV+^* (2ReA+aI)^+2Re^*jg^(£, t) 


If 0<a< Min ( | ReX (A) | ) 
1 ■ 


S3BP 


(23) 


Then 

-I2ReA+aI]-Q-Q’^>0 (A. 7) 

/. 0(z)— oV-W(z,t) (A. 8) 

Where 


W(^,t)"£*Q£-2Re£*^(£,t) (A. 9) 

Since £.(£,t) contains no terms linear in W(£,t)>0 provided 
I |z| I is sufficiently small. 

Thus 

V(z)^-aV for 1 1^| | sufficiently small 


Hence 

V(z) 6e"“*^V(0) (A. 10) 

2 I 2 

Thus if V(0)*-|||^|| is sufficiently small, V(z)«| | remains small and 

tends to zero as t-»-CD, hence (A. 2) is Liapunov asymptotically stable. Since 
JC|*T^, stability of ^ implies stability of x, hence (A.l) is L.A.S. 

Case (ii) A-defective, in this case it is not possible to diagonalize A. 

However, there exists a nonsingular matrix T such that A can be reduced to 
Jordan Canonical form 


t‘^at- 





(A. 11) 


(24) 


Where Che J are Jordan blocks such as 

o. 




0 



0 



(A. 12) 


To simplify the proof | consider the case 



Let x“T^ 
Then 


T"^AT« 


Xll 

0 X, 




dt 




0 



*+&(£♦ t) 



<A.13) 




I (A. 14) 





Let 


Where 


Where 


Since 


£(z,t)-T“^f(T*,t) 


I 1 8 (2 t ) 

Lim — * — “ — ^ »0 uniformly In t 

IklKo I UN 


y (A.u) 
I cont’d 




V(z)»z*t?z>0 


(A. 15) 


(ReX^)' 


T 

«P >0 


(A. 16) 


V«z*Pz+z*pz 


'-^*Q^+2Re (^*P£.(z , t) 


(A. 17) 
(A. 18) 



0 “2ReAjj_2 


(A. 19) 


ReX^<0 Vi 


T 

Q»Q >0 


(A. 20) 


Thus 


(26) 


Where 


W (£, t) "Z* (Q-al. ) z."2Re£*Pjg^(^, t ) 


(A. 22) 


If 

|ReX I 

a< y — » then (Q-al) > 0 (A. 23) 

(ReXj^)>l . 

Since S.(^,t) contains no terms linear In W(^,t), it is positive pro- 
vided I I is sufficiently small. 

V^-aV for I I^I I sufficiently small (A. 24) 

Applying the arguments of Case (l)» we see that the trivial solution of (A.l) 
is Liapunov Asymptotically stable. 

The technique developed above can easily be extended to cover the case 
of multiple repeated roots or higher order Jordan blocks. 


Critical Cases 

It will be observed that the techniques used to prove the stability of 
(A.l) break down if ReX^«0 for iG<l,k), i.e. if the matrix A has one or more 
zero eigenvalue, or one or more pairs of complex conjugate pure imaginary eigen- 
values. Such cases are called Critical Cases and will be treated in Theorem III. 

Theorem All 

If the matrix A in (A.l) has one or more eigenvalues with positive real 
part, then the trivial solution of (A.l) is Liapunov unstable for sufficiently 


small initial data. 


(27) 


Proof 

Caie (1) A nondeflectlve, in this case there exists a nonsingular matrix 
T which diagonalizes A. 

i.e., where 


Let x-T_z 
Then 


Let 


Where 


T”^AT-A 

ReX^>0 ie(l,k) 
ReAj SO j e (k+1) 


^ -Az+£(z,t) 

z ( 0 )m -t"^c 

i(z,t)-T"^f(z,t) 

Lim «o uniformly in t 

I I z| ho 1 1 z| i 

\h ° 1 



(A. 25) 


(A. 26) 


(A. 27) 


(A. 28) 


0 


(2BJ 


V»z*Pz+z*Pz 


V«r*>Qz+2Re (**Pfi.(£, t) ) 


(A. 29) 
(A. 30) 


Where 


2ReAk 


-2ReA 


N-k 


T 

•Q^SO 


(A. 31) 


Hence 


V-aV+W(z,t) 


(A. 32) 


Where 


W ( z , t ) [ Q-aP ]^+2Re^*P£(^, t ) 


(A. 33) 


If 0<a<Min X.(A) 
isi£;k ^ 


(A. 34) 


Then 


(Q-aP) is positive definite 


(A. 35) 


Since £.(^, t) contains no term linear in for 
W(^,t) is positive. 

Hence 

V^aV 


^ , sufficiently small, 


(A. 36) 


V(z)^e“*^V(0) 


(A. 37) 


Since V(z) is sign indefinite, there exists a set 


nj^;V(z)^0 dn^:V(z)=0 


(A. 38) 



(29) 


Define 


^ 2 * I 


Let 

^2! 0^002 (A. 39) 

If 

z(0)e«2 » V(0)>0 (A. 40) 

Then V(z)Se®^V(0)>0 (A. 41) 

The trajectory cannot exit through 30 since V(0)>0 and V(z) 

is increasing; therefore it must exit through the boundary ||^||”^< Hence 
given any 0<6<A there exists no e>0, such that if !!^(0 )||^e, £(0)£02* 
!l£(t)|l^fi ^or V t>0 the trivial solution of (A.l) is unstable in the 

sense of Liapunov. 

Case (ii) A defective, in this case A cannot be diagonalized, however there 
exists a nonsingular matrix T which reduces A to Jordan Canonical form 
i.e. 

k 

, I -N (A. 42) 

1«1 ^ 


T’^AT- 


al 


«2 


'J 


otk. 


Where 


(30) 


Let 


Case 


Where 


ai 


X,1 


Xll 


IJ 


<"Tz in (A.l) 


dz 

~«j£+£(z,t) 

£(0)-^-T“^£ 

J-T"^AT " 

£(z.t)-T"^f(Tz,t) 

Lim -0 uniformly in t 

Main ||z|| 

To simplify the presentation we shall consider three typical cases 


Xll 

0 X, 




ReXj^>0 , ReXj^O jG(3,N) 


(A. 43) 


(A. 44) 


(A. 45) 


(31) 


Let 


V(j8)»*^*P£ 


Where 


Then 



P- 


1 

(Re 



V-z^*Pz+js*pi 
■£*Qz+2Re£*P£(j^, t) 


Where 


Q-Q 


2ReXj^ 1 
1 2 


ReX, 


S 0 


-2ReA 


N-k 


equation (A. 49) may be rewritten 


Where 


V-aV+W(z,t) ; a>0 


W (£, t ) «£* (Q-aP ) +2Re^*P£^U , t ) 


If 

ReX- 

0<a< — 

l+(ReX^) 


(A. 46) 


(A. 47) 


(A. 48) 


(A. 49) 


(A. 50) 


(A, 51) 


(A. 52) 


(A. 53) 


(32) 


then 

<Q-aP) in positive definite (A. 54) 

Since £,(£,t) contains no terms linear in W(£,t) Is positive if 
||z||^A, sufficiently small. 


V^aV for 1 1 Z.I I ^ 4 


(A.55) 


Since V Is sign Indefinite, there exists a set 


Define 


From (A.55) 


VSiO , V-0 on (A. 56) 

(A. 57) 

nj^ns^2 (A. 58) 

V(z) e“‘^V(0) (A. 59) 


If 

} 

z(0)€fi3 ,9V(0)>0 

Then V(z)>0 and monotone Increasing provided I 1^1 16^3 • 

The trajectory, g^, starting In with V(0)>0 cannot exit through 
the boundary since V“0 on the trajectory must therefore exit through 

the boundary 11 ^ 11 “ 5 » Hence, given any 6 , 0<6<A, there exists no e>0, 
such that ||,z(0)|| '^e Implies | |£(t)|| s:6 V t>0. The trivial solution of 
(A.l) Is therefore unstable In the sense of Liapunov. 


(33) 


Cai« (lib) 


Where 


Lv?t 


Where 


In this case 




'kfl 


(A. 60) 


0 


‘N-k-2. 


ReXl>0 le(l.k) 

Rfi (^l^^j^) ^0 

ReAj^O j£(k+3,N) 


V(s)»£*Pz 




> (A. 61) 




(A. 62) 



-1 



1 

"^N-k-2 


(A. 63) 


V-£*Qz+2Rez*P£,(z , t) 


(A. 64) 


Where 


(3A) 


2ReA, 


0 


•2ReAk+i -1 


-1 


-Re\ 


k+1 


-2ReA 


N-k-2 


The matrix Q is clearly positive semi-definite. 
Equation (A. 64) may be rewritten 


(A. 65) 


V»aV+W(z,t) ; ct>0 


(A. 66) 


Where 


If 


W (£, t) (Q-aP )£+2Re^*P£^(£, t ) 


(A. 67) 


0<a<Min X. (A. 68) 
l^iaik 

The (Q-aP) is positive definite and W(^,t) is positive for |j^j|sA, 
sufficiently small. The arguments of Case (li) apply here also and the trivial 
solution of (A.l) is unstable in the sense of Liapunov. 


Case (lii) 



J- 



X 


z 


0 0 1 

0 oJ 


(A. 69) 


(35) 


Where 


Let 


Where 


Then 


Where 


ReX^>0 ie(l,k) 

je(k-»-l,N-2) 


j. (A. 70) 


V(z)«z*Pz 


(A. 71) 




0 



(A. 72) 


V“£*Q^-f2Re^*jg^(£, t) 


(A. 73) 


2ReAj^ 



“2ReA 

e 


0 -1 

-1 0 


(A. 74) 


Equation (A. 73) may be rewritten 


V-aV+W(z,t) ; 


o>0 


(A.75) 


(36) 


Where 

W (_z , t)^z* (Q-aP ) z+2Re^*P£(£, t ) (A . 7 6 ) 



0<a< Min X 



a 


y (A. 77) 

J 


Then (Q-aP) Is positive definite and W(£,t) is positive for 
sufficiently small. The arguments of Case (li) apply here also and the trivial 
solution of (A.l) is unstable in the sense of Liapunov. 

The techniques developed above are easily extended to the case of multiple 
repeated roots and higher order Jordan forms. 

It should be noted in passing that unlike Theorem I, Theorem II does not 
break down in the case where one or more eigenvalues have a zero real part. 


Critical Cases in the Liapunov-Poincare Theory 

As already pointed out, if the matrix A has any eigenvalues with zero real 
part, stability cannot in general be inferred from the stability of the linearized 
equations. In the case of the attitude stability of satellites it will be shown 
that due to the special form of the equations of motion, stability of the full 
perturbational equations can still be Inferred from the linearized or variational 
equations. 


Theorem AIII 

The perturbational equations governing the attitude stability of spinning 
satellites take the special form: 


( 37 ) 


3T = 


dx, 

“Ax2+|.2(Xj^.X2»t) 


Ai\ 

2L“I I 






uniformly in t 


S- (A. 78) 


If the matrix A is a stability matrix, then the trivial solution of (A. 78) 
is Liapunov stable for sufficiently small initial data. Furthermore, the states 
x^j^ and X 2 have the following properties 


Lim I Ix^l 1= 0 
t->oo 

Lim |l,x. ||«Y — constant 
t->-oo 


r (A. 79) 




Proof 

Let 


Where 




'v 


-2 




AT=A 


y (A. 80) 




ReX^<0 V i 



(38) 


We shall only discuss the case where A is nondefective, the case for 
A defective Is handled in a similar manner. 

Using (A. 80) equation (A. 78) becomes: 


where 



% (ll *^2 ’ ^ ''-1 ^-1 ‘ '^^2 ’ ^ 

£-2 ( *-2 ’ ^ » '^^2 ’ ^ 


''%(Zj^»0»t)=0 

I |£i(£i.£2‘*^^I I 

Llm =0 uniformly in t 

tiiaih Ikjll 




L (A. 81) 


Let 


V(z)=V^(z^)+V2(z2> 


(A. 82) 



(39) 


Where 



V2<l2>‘^2Vl 1^211^ 

j 

^ (A. 83) 

Hence 


V(£)-| l%l l^+l U 2 I |z| 1^ 

(A*84) 


^ 2 “— 2 ’ 

(A. 85) 

(A. 85) 

may be rewritten as 



V * ot'^O 

(A. 86) 

Where 


^2 (^1 »^2 * (Q“«I)^2"^*^®-2^2 ^“1 »~2 ’ ^ 

(A. 87) 

Where 


q=-2ReA>0 

(A. 88) 

If 


0<ot< Min |Re^j| 

(A. 89) 

Then 

. 



(Q-al) is positive definite 

(A, 90) 

Since 




/ 2 0 

.*. if 1 1^1 1” ^/ 1 1 +11^211 sufficiently small 

(A. 91) 


(AO) 


Then 

Where 

Thus by taking 

• • 

Hence 
Hence If 


From (A. 81) 
Since 

jniformly in t 


1 2 Re^ 5 B .2 

K 2 (A )*>0 as A->-0 


^ sufficiently small 


^ U 2 I 


V 2 ^-aV 2 




> (A. 92) 


(A. 93) 
(A. 94) 


V2(z 2> ^e”°“^V2(0) 


(A. 95) 


I U(0) I 1 ^ e<A 

^ 2 (^ 2 ^ ^ 

I \^2^t) 1 1 


(A. 96) 
(A. 97) 
(A. 98) 


.t 






(A. 99) 


^,^,t)= 0 and satisfies a nonlinearity condition; 


ll%(Zl»Z 2 ’''^ll ^Ki(A)V 2 (z 2 ) 

for ll^ll^A and Kj^(A)n.O(l) 




> (A. 100) 





Thus 




V2<Z2(t))dt 


(A. 101) 


Using (A. 97) 

I l%(t) I U I lli(O) I |+Ki<A) ” (l-e"“‘')e^ (A. 102) 

Using (A. 96) 


llz^Ct)!! ^K2(A)e 

(A. 103) 

V(z(t))-| |z(t)| lK 2 (A)+e"“*^]e^ 

(A. 104) 

1 |z(t)| 1 £;K3(A)e , V t>0 

(A. 105) 


Thus If we choose 


e ^ e*= 


6 

K3<A) 


(A. 106) 


Then given any 6, 0<6<A , if ] |^(0) [ | ^ e* then | \^{t) | j ^6 , V t>0 . 

Thus the trivial solution of the system (A. 81) Is Liapunov stable, and since 
(A. 81) is derived from (A. 78) by bounded linear transformations, system (A. 78) 
is also Liapunov stable. In addition, using equation (A. 98) we see that 

Llm ||z2(t)||=0 (A. 107) 

t->-oo 


Using (A. 80), (A, 107) implies that 



From (A. 99) 


Llm £j^(t)- 5 j^( 0 )+ £j^(Zj^,Z 2 ,T)dT 

t-^oo ' 0 


Since 


iij 


Ki(&) 2 

" e 

a 


» 


the integral (A. 109) converges, hence, using (A. 80) 


lim I !xi(t)i 1“Y “ a constant 

t->OD 


provided that the initial data is sufficiently small. 


Extension of Theorem AIII to Systems with Periodic Coefficients 

In the study of the attitude stability of dual spin satellites, the 
perturbational equations take the following form: 


I 

dF" “ 


dx^ 

^ ^“1 *~2 * ^ ^ 





(A. 109) 


(A. no) 


(A. Ill) 


A(t+T)=A(t) , V t 


( 43 ) 


llil(Xi»X2»*^^ll 

Llm 0 uniformly in t 

11*21 1 1*2! I 

For such systems, Theorem IV applies. 

Theorem AIV 

Given the system (A. Ill), if all the solutions of the equation 


dT i 2«»*‘=2 


are Liapunov asymptotically stable, then the trivial solution of (A, 111) is 
Liapunov stable for | ) sufficiently small. In addition the states 

and 3^2 have the following properties. 


Llm 11X9 11*0 

t*HX) 


Lim I |xj^| 1**Y ~ a constant 
t-^-ao 






J 


Proof 

Consider first the matrix equation 


^=A(t)X ; X(0)=I 

It is well known from Floquet theory that X(t) has the following form: 

X(t)=Q(t)e®^ 


Where B is a constant matrix 

and Q(t+T)=Q(t) , Q(0)=I Is a bounded periodic matrix. 


(A. Ill) 
Cont’d 

(A. 112) 
(A. 113) 

(A. 114) 

(A. 115) 
(A. 116) 



m 


The requirement that all solutions of equation (A. 112) be Liapunov 

I 

asymptotically stable Is equivalent to the requirement that the matrix B 
be a stability matrix, l.e, 

ReXj^(B)<0 , V 1 

The matrix Q(t) in (A. 115) satisfies the differential equation 


dCt 

dt 


+QB-A(t)Q 


Consider now the Liapunov transformation 


X2*Q(t)u2 


Substituting into (A. Ill) 


du 


^ U2+Q(t) — =A(t)Q(t\u2+f2(Xj^,Q(t)u2,t) 


dUrt . 

^ «Q"^(t)^ I A(t)Q(t)- 


A(t)Q(t)~ ^ u 




Using equation (A. 118) 


du. 




where 


.-1 


]l 2 »J ^2 * ^ *^^^-2 ^-1 ’ ^ ^ *^^-2 ’ 


Let 








- 1"~1 


(A. 117) 


(A. 118) 


(A. 119) 


(A. 120) 
(A. 121) 


(A. 122) 


(A. 123) 


(A5) 


Then (system (A. Ill) becomes 







■Bu2+h2<Uj^»U2'*^^ 



h^(u^,0,t)so ,l«l,2 


Ltm 

I IH2I 





“0 uniformly in 


t 


y (A.12A) 




The system (A. 124) has exactly the same structure as system (A. 78), hence 
by Theorem III, the trivial solution of (A. 124) is Liapunov stable for sufficiently 
small initial data and in addition Uj^ and have the following properties; 


i) Lim 

t-»-OD 

ii) Lim ] |u. (t) j 1 «y ““ a constant 
t->oD 


I (A. 125) 


Using (A. 119) and (A, 123) it therefore follows that system (A. Ill) is 
Liapunov stable and Xj^ and have the following properties; 


a) Lim ||x„(t)||*0 

t-^OD 

b) Lim I |x, (t) 1 |“Y “ a constant 

t->OD 


> (A. 126) 




Thus establishing Theorem IV. 


(46) 


AP PENDIX B 

JuBClElcatlon for the Method o£ Slowly VarylnR Parameters 

The linearized equations of motion of dual-spin sotelllteB with damping 
in both rotor and platform can be written in the following standard forms 



0 

1 




dx . . 
+A2(t)x 


] 


(B.l) 


^ ^ +D(t)v (B.2) 

Where v Is a two «;ector, is a four vector, 

D(t) are bounded matrices and s>0 Is a small parameter. 

To reduce (A.l) and (A. 2) to more convenient form, we Introduce the 
following transformation 


V" 


cos Xt sin Xt 
sin Xt -cos Xt 


(B.3) 


Equation (A.l) becomes 


Where 


d_a 

dx 


dx 


eA3<T) +eA^^(T)x 


A2(t)*T(Xt)Aj^(t) 

A^(t)“T(Xt)A2(t) 


T(Xt)= 


cos Xt 
oln Xt 


sin Xt 
- cos Xt 


(B.A) 




I (B.5) 




(47) 


Equation (B,2) bocomas 
where 

Let us now Introduce 


Substituting Into (B.6) 

Substituting (B.6) Into (B 
Where 

F(t). 

H(t) 


dx 


da 


d7 


(B.6) 


Cj^(t)-C(t)T(At) 
Dj^(t)-C(t) ~ +D(t)T(Xt) 


a second transformation 


> (B.7) 




x;»£+Gj^(t)^ 


I e “ Dj(5)d5 


-00 






(B.8) 


dz da 


(B.9) 


,A), using (B.8) and solving for 


dja 

dt 


da 

^ -eF(T)a+eHj^(T)£ 


(B.IO) 


.-1 


[I-eA3ej^] "IA3D^+(A3 Bq+A^)Gj^] 


.- 1 . 


[I-cAjC^l ‘UaBo+A^l 




Substituting (B.IO) Into (B.9) we have 


(B.12) 


Where 




*®02.+ Bj^(T)£feH2(t)a 


B^(T)»tCj(t)-G^(T)]Hj^(T) 

H2(t)*[C^(t)-Cj^(t)]P(t) 


y (B.13) 




For e sufficiently small, the matrices F(t), H 2 (t) and Bj^(t) 

are bounded 
if we write 

F(t)^Fq+Fj^(t) (B.14) 


where 


J 

F " Llm “ I F(x)d 
” T->co ^ h 


(B,15) 


Theorem BI 

Given the system of equations 

da 

-e[FQ+Fj^(T)]a+eHj^(T)£ 

dz 

~ «[Bq+ B^(T)l£feH2(T)a 
If 1) Fp and B^ are stability matrices 


■N 


r (B,16) 


11) 


Fj^(t)dT Is a bounded matrix 


(49) 


Then for e sufficiently smtill, the trivial solution of (B.16) Is ilapunov 
asymptotically stable. 


Proof 


Since Fjj and ate stability matrices there exist symmetric positive 


definite matrices P^ and such that 


Consider the function 




%Wo-h 


j^a+£^P 2£-' ( T ) £ 




I 


(B.17) 


(B.18) 


Where 


J 


Q(t)-F2(t)P^+Pj^F2(t) 

F2(t)“ I F^(t)dT 

Q(t) is a bounded matrix, since P^ and bounded matrices 

. »T T • »T T • I T* *T T 
V-a P^a+a P^a+z P 2 Z+Z la Qa+a Qa+a Qa} s 

Using equation (B.16) and (B.19) 


L (B.19) 


(B.20.) 


V-ea%Fj-»-Fj)Pj^+Pj^ (f'^(H‘QF) a 


4-2ea'^QHj^]+£'^ [ (Bp-4-EBj)P2+P2 <Bq+gBj^) ] z+2ea^H2P2Z 


A 


r(B.2i) 




(50) 


l^ere 

F+Fq+Fj^ 

Using (li.l7) V becomes: 


(B.22) 


V»-2ea'^a-£^^+e/(Bjp2+P2Bj^)£+2e A(t)£-e^ [£^(F'^(HQF)a+2a^^ (B, 23) 

Where 

R(t)«[P3^Hj^(t)+H2<t)P2J (B,2A) 

Equation (B.23) may be rewritten: 

V=-ea^a-/ [ I-e (B^p 2 +P 2 B]_) -eR^,R)£ 

- e (£-Rz ) ^ ( a,”R;3 ) 

-e^ (F"Q+QF) a+2a^QH^£] (B . 25) 

Since Bj^,P 2 >Q, R etc. are bounded matrices, for e sufficiently small, the 
sign of V is that of the first three terms. 

.*. V<0 for e sufficiently small (B.26) 


Similarly, the sign of V, (B.18) is that of the first two terms for e suf- 
ficiently small. Hence, for e sufficiently small 

V>0 , V<0 V is a Liapunov function (B.27) 

Hence, the trivial solution of equation (B.16) is Liapunov asymptotically stable. 

Using (B.8), stability of £ and a^ implies stability of jc and 
and hence of x and v. Thus, under the hypothesis of Theorem (BI) , the 


(51) 


trivial oolutlon of equations (B.l) and (B.2) is Liapunov asymptotically stable. 

Given that B^ is a stability matrix, the requirements for stability are 
that e be sufficiently small and that the time average of the matrix F(t) 
should be a stability matrix. 

Non? 


F(T)-[I-eA3(T)Cj^(T)j"^[A3(T)BQ-fA^(T))Gj^(T)+A3(T)Dj^(T)] 

■ (A3(t)Bq+A^(t))Gj^(t)+A3(t)D^(t) 

+cII-eA3Cj^]’^A3C^( (A3 Bq+A^)Gj^+A3Dj^ (B.28) 


Hence 


F *F +eF 
'^0 *^00 ^®01 




Where 


,00 


Foo=Llm ^ 

T-+00 


[ (A3BQ+A^)Gj^+A3D]d 


y (B.29) 


,a> 

f L [I-^A C r^AjC^tCAjE +A )G +A D ]dr 
T'+oo ^ 0 


The requirement for stability is that Fq be a stability matrix, for 
e sufficiently small this requirement will be satisfied if the matrix F^^ 
in (B.29) is a stability matrix. 

In terms of the matrices A^, A2, B^, C and D, this condition becomes: 


(52) 



wLlm 

T-^oo 


.T 


r'^ Bq(t-0 
^0 



T(XOde 


+Aj^[C ^ +D]T(Xt) 


(B.30) 


should be a stability matrix. 

Rather than carry out the operations In (B.30) In one step, may 

be evaluated In the following manner. 

i) Make the transformation (B.3) 

11) Compute 2io('‘‘) "" "steady state" response of equation (B,6) 
regarding £(t) as a constant vector 
ill) Substitute the "steady state" response ^(x) Into equation (B.4) 
iv) Time average equation (B.4) regarding _a(t) as a "constant" vector 
The procedure yields the equation 


dt 


=^00^ 


(B.31) 


Where is exactly the expression In equation (B.30). 

It will be noted that this procedure is exactly what was done heuristically 
In Section 2 of this note. 
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APPENDIX C 

Equations of Motion for a Dual-Spin Satellite 

The dual-spin satellite consists of two rigid bodies having a common 
axis of rotation. 

Let the axial moment of Inertia of the rotor be J. 

Let the total axial moment of Inertia of the satellite be C (Including platform, 
rotor and dampers) 

Let the total equatorial moment of Inertia of the satellite be A (including 
platform, rotor and dampers) 

To simulate the effect of Internal damping, the model will include torsional 
dampers. 

Let i|) be angle (about the common axis) between the body fixed axis in the 
rotor and platform. 

Equation of Motion 
Rotor 

]=T^-Tg (C.l) 

Where l!. is the polar moment of inertia of the damper on the rotor 
b 

0 is the rotation angle of the torsional damper on the rotor 

T^ the torque nf the despin motor 

T^ the frictional torque of the rotor bearings. 

B 



(54) 


Rotor/Platform 


-03j^sin i|j )-l|^<j)W2*0 




AWj^-f [ (C-A)w2+J<I']w2+I^[^cos ij; -G(w2+'j^)sin tp ]+I^$*0 


> (C.2) 


Aw2“((C“A)w2+J'i']Wj^+I^[0sin i|) +0(w2+ti')cos ip 


• . • . 




Where is the polar moment of inertia of the torsional damper on the platform 

<(i is the rotation angle of that damper * 


Damper s 


I^e+K^e+K20=-l|j[(I)j^cos t|j +o)2sln ip +o)2^cos ip -ojj^ij/sin ^ ] 

I^6+K^J+K2 


Steady State Solution 


a>^=fi , ip‘=0 , iji=at 
0=(j)=6=(|»=Wj^=W2^ 0 


Perturbed Motion 


Let 


Q 


ti)2=J^+a^ ; T=ot ; r= — 


i|;= \p2=^-^=o(l+z) 


0) w 

— =v- ; — =v, T_-T =gac 

a l a 2 B M 


oMuaBisiiiaisaiaaysaKSii 


y (C.3) 








(C.4) 


(C.5) 
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Perturbatlonal Equations 


[3^ + df] [v2Cos(x+n)-ViSln<T+n)] + fc -» 


dx ^ 


^ C ^ [v2Cos(t+n)-v^8in(T+n)] ^ -W I V2 ^ “0 

^ + |^x+ ^ +m' 1^^ cos(x+n)“ ^ <l+r+g+c)sln(x+n)j “0 

^ _ j^x+ j Vj^+M(r+C) ^ +\i' sin(x+n)+ — (l+r+5+5 )cos(x+n) 


d e de 

TT dx 

dx 


26+ r +V2(1+?)) cos(T+n)+ -v^(l+!;)^ sin(T+n)J “0 


4 +*=1 f ^ -0 

dx 




(56) 


Let 



0 

ISj 


Using (C.8) and (C.9), equations (C.6) can be rewritten as 5 


=A2E+M[B^(T)+B2(i,£,T)] ~ +i2(i»E»T) 


(C.8) 


(C.9) 


(C.IO) 


(C.ll) 


Where 


A^, A2 are constant matrices 
Bj^(t+2tt)=B^(t) 

Lltn — =0 uniformly in t 

i lil l+l |£| N IUIIh-IIbII 




y (C.12) 


4(1,0,t)=0 

i2(l,0,T)=0 




Bnaaiiii 






Llm = -0 

I IeI I-o I W I 


I Ii2(i»£»^> 1 1 

Iliil+llElho llelkIWl 

Solving (C,ll) for d£/dT, we havej 

dp 

“A3(x)£+f3(i,£,T) 


Where 

A3(t+7t)“A3(t) 

+pn-p(B^(T)+B2(i,£,T) r\(i.£, t)A3£ 


Where 


i3(ii0,T)“j0 


and Llm 


li‘j(i.£.'>^) 1 1 


111 Kl l£l llill+ll£ll 


=0 uniformly in t 


Let 


where T ^AT= 


0 0 

0 

0 -y 


y>0 




V (C.12) 
I cont’d 




(C.13) 


y (C.14) 




> (C,15) 




(C.16) 
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Equations (C.IO) and (C.13) become 


dt 


d£ 


-YJ 




^ “A3(T)£ff5(X,£,T> 


> (C.17) 


Let 




Equations (C.17) can now ba written in the form; 
dx^ 




(C.18) 


^ »A(T)x2+f2(Xi.X2‘'^^ 


x» 



x(0)=c 


A(t+it)®A(t) 
f^(x^,0,T)«0 i=l,2 

I I 

Lim ; — =0 uniformly in t 

llxJh*^ Mi? I 


(C.19) 


(59) 


» 


Whore 



(C.20) 


Theorem IV of Appendix A is applicable to (C.19) and stability is 
guaranteed, for | | c| | sufficiently small, provided all solutions of 


dj«2 

^“A(t)x2 (C.21) 


are Liapunov asymptotically stable. 

r 

Using (C.20), all solutions of (G.21) will be asymptotically stable 
provided all solutions of 

d£ , V 

^ “A3(t)£ (C.22) 

are asymptotically stable. 

Since equation (C.22) is the linear part of equation (C.13), the condition 
for stability can be expressed in terms of the variational equations obtained 
by linearizing the perturbatlonal equations (C.6), 

Thus, the conditions for stability of a dual-spin satellite are: 

1) The initial perturbations shall be sufficiently small 

2) All solutions of the following set of differential equations shall 
be asymptotically stable 


< 60 ) 



It will be noted that the present analysis rigorously justifies the normal 
engineering analysis, in which one examines only the stability of the linearized 
equations and Ignores completely the subtleties of the stability of the per- 
turbational equations. 


i 


(6X) 
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